Bimodal Counting Statistics in Single Electron Tunneling through a Quantum Dot 
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We explore the full counting statistics of single electron tunneling through a quantum dot using a 
quantum point contact as non-invasive high bandwidth charge detector. The distribution of counted 
tunneling events is measured as a function of gate and source-drain-voltage for several consecutive 
electron numbers on the quantum dot. For certain configurations we observe super- Poissonian 
statistics for bias voltages at which excited states become accessible. The associated counting 
distributions interestingly show a bimodal characteristic. Analyzing the time dependence of the 
number of electron counts we relate this to a slow switching between different electron configurations 
on the quantum dot. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.Kv 



The dynamics of electron transport through a quan- 
tum dot cannot be accessed by measurements of the av- 
erage (DC) current alone. Additional information was 
successfully gained from measurements of the shot noise 
P, 0, Q and recently even of the 3rd moment [4] of the 
current through the quantum dot. But it is hard to see 
how such measurements could be extended to even higher 
moments. 

Recently an alternative approach was pointed out for 
semiconductor quantum dots: Beside a direct measure- 
ment of the resonant tunneling current through the dot 
one can also measure the charge on the quantum dot us- 
ing a nearby quantum point contact as a non-invasive and 
highly sensitive detector 0, @, 0, Q . For sufficiently low 
tunneling rates and high detector bandwidth this allows 
to resolve individual tunneling events onto and off the 
dot [3,113 and to measure the full counting statistics for 
tunneling through a quantum dot ground-state [ll|, [l2[ . 

We have implemented a measurement of single electron 
counting with large bandwidth. This enables us to mea- 
sure the full counting statistics for single electron trans- 
port through a quantum dot as a function of both bias 
and gate voltage for a series of consecutive electron num- 
bers on the dot. We observe well understood counting 
distributions for tunneling through the ground-states of 
the quantum dot. But interestingly our analysis reveals 
also the occurrence of bimodal counting distributions for 
certain numbers of electrons on the dot and sufficient bias 
voltage. We relate this to a slow switching between two 
different quantum dot configurations that have distinct 
tunnel couplings to the leads. 

Our device is based on a GaAs/AlGaAs heterostruc- 
ture containing a two-dimensional electron system 
(2DES) 34 nm below the surface. The electron density is 
n = 4.59 • 10 15 m~ 2 , the mobility is /i = 64.3 m 2 /Vs. We 
have used an atomic force microscope (AFM) to define 
the quantum dot (QD) and the quantum point contact 
(QPC) structure by local anodic oxidation (LAO) on the 
surface pjj 03 ; the 2DES below the oxidized surface is 
depleted and insulating areas can be written. 

An AFM image of our device is presented in Fig. [IJl 
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FIG. 1: Operating principle of the device containing a QD and 
a QPC. (a) AFM image of the device and gate configuration, 
(b) Time segment of the QPC signal, tunneling times Ti n and 
Tout are marked, (c) Distribution of tunneling times extracted 
from 100 time segments. Each time segment spans 5 ms. The 
distribution for n n is shifted by a factor of five for clarity, 
(d) Distribution of tunneling events for same measurement 
compared to a theoretical model calculation. In total 446 
electrons passed the QD during the measurement leading to 
this distribution. 



The bright walls depict the insulating lines written by 
the AFM. The QPC (left area) is separated from the QD 
structure (right area) by an insulating line. The QPC 
can be tuned using the in-plane gate G3. The QD is 
coupled to source and drain via two tunneling barriers, 
which can be separately controlled with gates Gl and G2. 
These gates are also used to set the number of electrons 
in the QD. We use two electrically separated circuits 
for simultaneous conductance measurements through the 
QPC and the QD. The experimental setup allows us to 
resolve tunneling times r as small as 30 ills. In Fig.PJ) the 
time resolved current through the QPC is shown. Dis- 
tinct changes in the current occur for every change of 
the number of electrons on the nearby QD. For exam- 
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pie the current drops down whenever an electron enters 
the dot and thus the number of electron changes from 
N to N+l. Due to the sufficient large bias voltage of 
0.2 mV electrons are entering the QD solely from the 
source and are leaving the QD to the drain. In the time 
segment shown here five traversing electrons can be iden- 
tified. Additionally the tunneling times T in and r ou t can 
be extracted from the QPC signal where T in is the time 
till an additional electron hops from source onto the dot 
and T out is the time that an electron needs to leave the 
dot into the drain. 

The statistical distributions of the tunneling times r- in 
and Tout are shown in Fig. [it. They follow an exponential 
decay, where the exponent is given by the tunneling rates 
Tin and r ou t respectively. We find Ti n = 1.7 kHz and 
r ou t = 1-9 kHz. Alternatively one can also determine 
these rates directly from the mean value of the tunneling 
times using the relation r in ( out ) = (r in ( out ) We find 
a good agreement between both methods if the rates are 
small compared to the bandwidth of our measurement. 
For fast rates we will miss a number of short transitions, 
which affects the mean of the tunneling times [Hj]. Here 
the analysis of the slope of the distribution yields superior 
results. 

In this paper we will additionally apply a further 
method to extract the statistical properties of the sin- 
gle electron transport. The measured long time trace is 
divided into a large number of short segments as depicted 
in Fig. [TJd and the number of transitions n from N+l to N 
electrons is counted for each segment. We now determine 
the statistical distribution of these counts n as depicted 
in Fig[T]i. The experimental results (bars) compare well 
to the theoretical model (line) calculated without free pa- 
rameters from the mean of r- in and r ou t [HI Ha]- From 
the distribution one can extract not only the mean value 
of n but higher moments like the Fano factor a which is 
given by the second moment of the obtained distribution 
divided by the mean value (n): 
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Now we will study the statistics in detail for nonlinear 
transport through the QD. In Fig. [2^ the rate of elec- 
trons per second traversing the QD as deduced from (n) 
is shown as a function of gate voltage and bias. Each 
data-point represents a full measurement of the distribu- 
tion of counting events and the rate is determined from 
the mean value (n). We clearly observe the so-called 
Coulomb diamonds well known from conventional trans- 
port experiments. Clear Coulomb blockade regions are 
found as well as discrete regions of finite current due to 
single electron transport through the ground state and 
for large bias also through excited states of the QD. The 
QPC detector can also be used to determine the mean 
charge of the QD. [17] For this we analyzed the DC cur- 
rent through the QPC and extracted the mean charge in- 
formation. The changes in the mean charge of the QD are 
shown as lines in Fig.[2|). The DC charge detection corn- 
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FIG. 2: (a) Counting rate of electrons passing the QD as a 
function of gate voltage and source-drain voltage. The lines 
mark changes of the average QD occupation which are de- 
duced from the DC-average of the detector signal, (b) Fano 
factor a extracted from the full counting statistics for the 
same range. For the white areas above -130 mV at least one 
of the tunneling rates is getting to high to extract a reliable 
counting statistic. 



pares well with the results of the real-time measurement. 
When the mean charge of the dot changes a distinct step 
in the counting rate can be seen. 

In Fig. [2)3 the Fano factor a as determined from Eq. 
[T]is shown for the same measurement. It turns out that 
a is between 0.5 and 1 for most of the QD configura- 
tions as one would expect for single electron transport 
through the ground state of a QD [3]. But for a spe- 
cial area marked by the red triangle in Fig. [2] (region A), 
super-Poissonian noise is observed. In the marked range 
not only the ground-state but also excited states take 
part in transport. While the ground- state transport be- 
low region A shows sub-Poissonian characteristic, a rises 
dramatically when an excited state enters the transport 
window. A similar behavior can also be observed for the 
opposite transport direction at the next higher electron 
number (Fig. O region B). Below region B the ground- 
state transport is strongly suppressed but when the ex- 
cited state takes part in the transport the current rises 
and an a > 1 is observed. 

We now analyze the counting statistics for region A 
marked by the triangle. For finite bias but only the 
ground-state inside the bias-window (directly below re- 
gion A), we find no special behavior. Applying the 
same procedure as described for Fig. []J we obtain tun- 
neling rates of r out = 3.4 kHz and Ti n = 2.8 kHz for 
V s d = 1-48 mV and Vq = —143 mV . But as we cross 
into region A, i.e. as an excited state enters the bias 
window, we find a completely changed counting statistic 



3 



with a bimodal distribution. A characteristic distribu- 
tion is shown in Fig. [3^. In the distribution two clearly 
distinct maxima can be identified instead of the single 
peak structure of the usual distribution function. Inter- 
estingly one of the maxima emerges at a lower number of 
counts than the mean (n) observed for tunneling through 
the ground-state before an excited state comes into play. 
How does the distribution relate to the rate and the Fano 
factor displayed in Fig. [2]? Firstly the bimodal structure 
leads to a significantly higher width of the distribution 
that results in the strong increase of the Fano factor. In 
contrast the mean value does not show such a dramatic 
change (compare Fig. [2^i), whereas the mean value of the 
bimodal distribution itself has a fairly low probability to 
be observed. 

To study the origin of the bimodality we now analyze 
the evolution of the counted electron events with time. 
This is shown in Fig. [3)3. Depicted is the number of 
electrons entering the QD in a given time segment in 
chronological order. By this we can examine the tempo- 
ral evolution of all time segments contained in the distri- 
bution of Fig. [3^. As expected from the bimodal distri- 
bution there are two separated ranges where the major- 
ity of events are accumulated. These are marked by the 
two darker bands in the back of Fig. [3](a,b). Moreover 
a switching between these two ranges can be observed. 
Typically the system stays in each of this two ranges for a 
time r s of 50 to 250 milliseconds what is fairly long com- 
pared to the typical tunneling times (r in / out ) of about 0.2 
to 1 milliseconds. For the second region showing super- 
Poissonian noise (Fig. [2j region B) we found the same 
bimodal characteristic and also a comparable switching 
behavior. Therefore this bimodality seems to be a more 
general feature observable at different electron numbers. 

Our observations can be explained by the existence of 
two different QD configurations. Each configuration pro- 
vides a transport channel for electrons. This behavior 
leads to super-Poissonian noise, but differs from the sit- 
uation described in [H, Qjl for multilevel quantum dots. 

To acquire a deeper understanding of the ongoing pro- 
cesses we analyzed the tunneling times of the two trans- 
port channels in detail. For this we individually studied 
tunneling events that can be related to one of the two 
QD configurations. 

The resulting distributions of the tunneling times 
r in(out) are shown in Fig. [3]i separately for configuration 
I (1st and 3rd graph) and II (2nd and 4th graph). All 
sets of data were fitted by exponential decay to extract 
the tunneling rates. While the tunneling barriers were 
roughly symmetric for the ground state as deduced from 
similar values of Ti n and r out (see FigJT] and the values 
given for V s d = 1.48 mV), they get clearly asymmetric 
for the analyzed situation. 

When the first excited state enters the transport win- 
dow, the tunneling rate Ti n has risen rather evenly for 
both configuration I (6.0 kHz) and II (5.7 kHz). Both 
configurations seem to differ only slightly in the coupling 
to the source lead. For the drain lead the result is differ- 
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FIG. 3: (a) Typical bimodal distribution of electrons passing 
through the dot for the area marked by a triangle in Fig. [2] 
Transport channels are marked by a number and a different 
background (color online). The experimental results (bars) 
compare well to the basic model calculation (line), (b) Time 
evolution of the detected number of electrons per time seg- 
ment. A switching between the two transport channels occurs, 
(c) Simple scheme of the experimental findings, (d) Distribu- 
tion of the tunneling times r in ( out ) for both transport channels 
(dots). Short tunneling times (light gray) are suppressed due 
to limited bandwidth. The denoted tunneling rates r in ( out ) 
are extracted from the exponential fit (line). 



ent. r out for configuration I (3.6 kHz) is still comparable 
to the ground-state. In contrast configuration II shows a 
tunneling rate r out of only 1 kHz that is far less than the 
outgoing tunneling rate of the ground-state transport. 
When the QD is in configuration II tunneling off the dot 
is suppressed as the QD seems to couple less efficiently 
to the drain lead. 

We can model the strong bimodality of the counting 
distribution assuming two independent transport chan- 
nels. Starting from the full counting statistics theory [l2[ 
for single state transport through a QD we combine two 
distribution functions to model our experimental results. 
The probability to measure a number n of traversing elec- 
trons is given by: 

P(n) = q J\- s ^-^ + (l- q ) j\s 2 U)-n ( ^ (2) 
where Si and £2 are the generating functions for the sep- 
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arate transport channels and q is the probability to detect 
transport through channel I. The generating function of 
a single level QD has been calculated for unidirectional 
tunneling as observed at sufficient large bias [Hj]: 

^ = (r in + Tout - \j (r in - r out ) 2 + 4r in r out e— ^ (3) 

We use the tunneling rates extracted individually for 
the two QD configurations (Fig. [3]i) which leaves only q 
as remaining free parameter. The resulting distribution 
function shown in Fig. [HJa, dashed line) was received for 
q = 0.51. The good agreement of model and experimen- 
tal data further confirms the idea of a switching between 
two dot configurations where the dynamics of each con- 
figuration individually is described by a single pair of 
tunneling rates r in and r out . 

The most likely nature of the different configurations 
is the excitation of an electron into a non-equilibrium 
single particle state that couples only weakly to source 
and drain. If the ground-state transition is given by 
(N) -> (TV + 1) -> (N) we could relate this to config- 
uration I. Thus the change of configuration I to II would 
involve a process (N) -> (N + 1*) -> (N*) with (N) and 
(N + 1) the ground states for N and N+l electrons and 
(N*) an excited state. If the lowest lying free single par- 
ticle state is only weakly coupled, we can carry a single 
electron current in the cycle (N*) — > (N + 1*) — > (N*) 
(compare Fig.[3t) until the less probable (N + 1*) — > (N) 
transition reinstates configuration I (ground-state tran- 
sition). 



The two configurations can arise from two origins: spin 
or charge. We find that the bimodality only appears 
when excited states with a significant excitation energy of 
0.3 meV are accessible. We assume that the occurrence of 
two transport channels is caused by two different charge 
configurations. A change of the spin configuration should 
mainly change the tunneling rate into the dot due to a 
change of spin selection rule and should have only small 
effect on the outgoing rate. In contrast we observe for 
both configurations roughly the same tunneling rate for 
electrons entering the QD, while a significant difference 
occurs in the outgoing rate. The experimental results 
therefore favor a charge type effect. 

To conclude, we have shown full counting statistics of 
single electron tunneling through a quantum dot using a 
quantum point contact as non-invasive high bandwidth 
charge detector. We observe super-Poissonian noise for 
certain QD configurations where excited states take part 
in electron transport. For these configurations a clear 
bimodality of the electron counting distribution occurs. 
We analyzed the bimodal distribution in detail and found 
a slow switching behavior in the dot transport. We an- 
alyzed the tunneling times for both configurations and 
presented a model. The good agreement of model and 
experimental results confirms the presumption of two in- 
dependent transport channels. The switching between 
two independent transport channels can be explained by 
the existence of two different QD configurations. 

This work has been supported by BMBF in the frame- 
work of nanoQUIT. 
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